MGGPOD: A MONTE CARLO SUITE FOR GAMMA-RAY ASTRONOMY 
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ABSTRACT 



We have developed MGGPOD, a user-friendly suite 
of Monte Carlo codes built around the widely used 
GEANT (Version 3.21) package. The MGGPOD 
Monte Carlo suite and documentation are publicly 
available for download. MGGPOD is an ideal tool for 
supporting the various stages of gamma-ray astron- 
omy missions, ranging from the design, development, 
and performance prediction through calibration and 
response generation to data reduction. In particu- 
lar, MGGPOD is capable of simulating ab initio the 
physical processes relevant for the production of in- 
strumental backgrounds. These include the build-up 
and delayed decay of radioactive isotopes as well as 
the prompt de-excitation of excited nuclei, both of 
which give rise to a plethora of instrumental gamma- 
ray background lines in addition to continuum back- 
grounds. 

Key words: gamma-ray astronomy; Monte Carlo 
simulation; instrumentation. 



1. INTRODUCTION 



Intense and complex instrumental backgrounds, 
against which the much smaller signals from celestial 
sources have to be discerned, are a notorious prob- 
lem for low and intermediate energy gamma-ray as- 
tronomy (~50 keV - 10 MeV). Therefore a detailed 
qualitative and quantitative understanding of instru- 
mental line and continuum backgrounds is crucial for 
most stages of gamma-ray astronomy missions, rang- 
ing from the design and development of new instru- 
mentation through performance prediction to data 
reduction. A promising approach for obtaining quan- 
titative estimates of instrumental backg rounds is ab 
initio Monte Carlo simulation (see e.g. iDean et all 
12001 . 
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We have developed a suite of Monte Carlo pa c kages , 
name d MGGPOD l|Weidenspointner et all 1200.1 
2004), that supports this type of simulation. The 
MGGPOD Monte Carlo suite (version 1.0) and doc- 
umentation are publicly available for download from 
http://sigma-2.cesr.fr/spi/MGGP0D/. In this 
paper we provide an overview of the capabilities, 
functioning, and structure of the MGGPOD pack- 
age, and give examples of applications to past and 
present gamma-ray missions. 



THE MGGPOD MONTE CARLO 
SIMULATION SUITE 



The MGGPOD Monte Carlo suite allows ab ini- 
tio simulations of instrumental backgrounds - in- 
cluding the many gamma-ray lines - arising from 
interactions of the various radiation fields within 
the instrument and spacecraft materials. It is pos- 
sible to simulate both prompt instrumental back- 
grounds, such as energy losses of cosmic-ray par- 
ticles and their secondaries, as well as delayed in- 
strumental backgrounds, which are due to the de- 
cay of radioactive isotopes produced in nuclear in- 
teractions. Of course, MGGPOD can also be used 
to study the response of gamma-ray instruments to 
astrophysical and calibration sources. The MGG- 
POD suite is therefore an ideal Monte Carlo tool 
for gamma-ray astronomy. A detailed description of 
the physics of the M GGPO D Monte Carl o suite can 
be found in Weidenspointner et al. (2004). The doc- 
umentation available from the MGGPOD web site 
provides comprehensive practical advice for users. 
Fi rst applications of MGGPOD have b een presented 
in lWeidenspointner et all l|2003Ll2004[ l. 



2.1. Capabilities and Functionalities 

MGGPOD is a suite of five closely integrated Monte 
Carlo packages, namely MGEANT, GCALOR, 
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Figure 1. A flow chart illustrating the overall structure of the MGGPOD Monte Carlo simulation suite. The 
various simulation packages (shown in boxes) and input and output files (shown in ellipses and round- edged 
boxes) are explained in the text 



PROMPT, ORIHET, and DECAY. The MGG- 
POD suite resulted fro m a co mbinati on of the 
NASA/GSFC MGEANT l)Sturner et a l.. 200(3) and 
the U niversity of Southampton GGOD l|Dean et all 
2003) Monte Carlo codes, which we supplemented 
with the newly developed PROMPT package. All 
these packages are based on the widely used GEANT 
Detector Description and Simulation Tool (Version 
3.21) created and supported at CERN 1 , which is de- 
signed to simulate the passage of elementary particles 
through an experimental setup. 

In a nutshell, the capabilities and functions of the five 
packages that constitute the MGGPOD suite are as 
follows: 

• MGEANT is a multi-purpose simulation pack- 
age that was created to increase the versatil- 
ity of the GEANT simulation tool. A modular, 
"object oriented" approach was pursued, allow- 
ing for rapid prototyping of detector systems 
and easy generation of most of the radiation 
fields relevant to gamma-ray astronomy. Within 
the MGGPOD suite, MGEANT (i.e. GEANT) 
stores and transports all particles, and treats 
electromagnetic interactions from about 10 keV 
to a few TeV. MGEANT provides the option of 
using the GLECS package 2 to take into account 
the energy of bound electrons in Compton scat- 
terings. The MGEANT simulation package and 
a user manual are available at a NASA/GSFC 
web site 3 . 



• GCALOR (|Zeitnitz and Gabriel Il994() simu- 
lates hadronic interactions down to 1 MeV for 



1 see http://wwwinfo. cern.ch/asd/geant/ 

2 see http://nis-www.lanl.gov/~mkippen/actsim/glecs/ by 
R.M. K ippen 

3 see http://lheawww.gsfc.nasa.gov/docs/gamcosray/legr/ 
mgeant/ 'mgcant.html 



nucleons and charged pions and down to ther- 
mal energies (10 -5 eV) for neutrons. Equally 
important, this package 4 provides access to the 
energy deposits from all interactions as well as 
to isotope production anywhere in the simulated 
setup. 

• PROMPT simulates prompt photon emission 
associated with the de-excitation of excited nu- 
clei produced by neutron capture, inelastic neu- 
tron scattering, and spallation. 

• ORIHET, originally developed for the GGOD 
suite and improved for MGGPOD, calculates 
the build-up and decay of activity in any sys- 
tem for which the nuclide production rates are 
known. Hence ORIHET can be used to con- 
vert nuclide production rates, determined from 
simulations of cosmic-ray irradiation, to decay 
rates. These are required input for simulat- 
ing the radioactive decays giving rise to delayed 
background. 

• DECAY, again originally developed for GGOD 
and improved for our purposes, enables MGG- 
POD to simulate radioactive decays. 



2.2. Structure 

The overall structure of the MGGPOD package is 
illustrated in Fig. ^ Depending on the simulated 
radiation field or gamma-ray source distribution one 
or three steps, requiring two or three input files, are 
needed to obtain the resulting energy deposits in the 
detector system under study. In general, it is advis- 
able to simulate each component of the radiation en- 
vironment separately. MGGPOD distinguishes two 

4 see http://ws www.physik. uni-mainz. de/zeitnitz/gcalor/ 
gcalor.html 



classes of radiation fields. Class I comprises radia- 
tion fields for which only prompt energy deposits are 
of interest, such as celestial or laboratory gamma- 
ray sources or cosmic-ray electrons. Class II com- 
prises radiation fields for which in addition delayed 
energy deposits resulting from the activation of ra- 
dioactive isotopes need to be considered. Examples 
for Class II fields are cosmic-ray protons, or geomag- 
netically trapped protons. 

For both of these classes, the simulation of the 
prompt energy deposits requires two input files: a 
mass model, and a model of the simulated radiation 
field. The mass model is a detailed computer de- 
scription of the experimental setup under study. It 
specifies the geometrical structure of instrument and 
spacecraft, the atomic and/or isotopic composition 
of materials, and sets parameters that influence the 
transport of particles in different materials. Each 
component of the radiation environment (and anal- 
ogously for gamma-ray sources) to which the instru- 
ment is exposed is characterized by three quantities: 
the type of the incident particles, and their spectral 
and angular distributions. The prompt energy de- 
posits are written to an output event file; in case of 
a Class II radiation field there is an additionial out- 
put file in which all the nuclei produced in hadronic 
interactions are recorded. 

To simulate delayed energy deposits (Class II radi- 
ation field) two additonal steps need to be taken. 
These require as input the time history of the radia- 
tion field which is responsible for the activation, and 
the previously calculated isotope production rates. 
Based on this information first the activity of each 
isotope produced in each structural element of the 
mass model is determined. Then these activities are 
used to simulate the delayed energy deposits due to 
radioactive decays in the instrument. 

Combining prompt and delayed energy deposits from 
each component of the radiation environment and 
gamma-ray sources, it is possible to obtain the to- 
tal energy deposited in the system as a function of 
position and time. 



3. PRACTICAL CONSIDERATIONS 



This section addresses some of the practical consid- 
erations a potential user of MGGPOD might have. 
It is a synopsis of the documentation available from 
the MGGPOD web site. 

Installation of the complete MGGPOD package 
(software and data files) requires little more than 
100 MB of disk space. The documentation provides 
detailed installation instructions, including the in- 
stallation of the required CERNLIB and other li- 
braries which are necessary to build, but are not in- 
cluded in, the MGGPOD package. The MGGPOD 
software is written in the FORTRAN 77 and C pro- 
gramming languages; in addition there are a few C 
shell scripts. The source code is completely open 



and transparent, allowing the user to adapt, change, 
and improve the code. In fact, a few customiza- 
tions and adjustments of the code are inevitable 
when simulating different instruments. These un- 
avoidable changes are described in the documenta- 
tion. The format of both the event and activation 
output files is FITS. MGGPOD inherited the inter- 
active display capabilities of (M)GEANT, which are 
based on CERN's PAW++ package. This is a very 
convenient feature when creating a mass model, or 
when specifying the parameters of a radiation field 
("beam"). 

To facilitate the use of MGGPOD by novice users 
the release contains examples for all simulation steps 
outlined in Sec. 12.21 The examples include all input 
and output files, and instructions on how to build 
and run the necessary executables. 

As described in Sec. 12.11 different modules of MGG- 
POD simulate different physical processes which can 
be treated over different energy ranges. The MGG- 
POD input file allows the user to define low-energy 
cutoff values for the tracking of five different parti- 
cle types: photons, electrons, neutral and charged 
hadrons, and muons. We recommend to use as de- 
fault low-energy cutoff values 10 keV for photons, 
electrons, and muons, 10~ 2 eV for neutral hadrons, 
and 1 MeV for charged hadrons. Experienced users 
can lower these thresholds, in particular for photons. 

For many applications of MGGPOD the computa- 
tion time needed to complete a simulation is an 
important consideration. The speed of MGGPOD 
Monte Carlo simulations varies greatly with the type 
and energy of the incident particle (which affect e.g. 
the number of secondary particles that are produced 
and tracked) , with the physical volume and mass de- 
scribed in the mass model (but less with geometrical 
complexity) , and of course with the speed of the com- 
puter employed. As an example, for the TGRS and 
SPI mass models, using a 1.6 GHz CPU, typically 
about 3 x 10~ 3 CPU seconds and 2 x 10~ 2 CPU sec- 
onds are needed to simulate the interactions of a 
cosmic-ray proton, respectively. The difference in 
processing time is mainly due to the fact that LNTE- 
GRAL/SPI is a much heavier and complex instru- 
ment than WIND/TGRS. For both models the time 
needed to simulate an incident gamma-ray photon, or 
a radioactive decay, is about ten times shorter than 
the respective proton processing time. The compu- 
tation time for a simulation also depends on the de- 
sired statistics in the result; better statistics require 
the simulation of more particles and consequently 
more processing time. For SPI, about 10 7 photons, 
requiring a computation time of about 8 CPU hours, 
in a homogeneous beam covering the 19 Ge detec- 
tors are sufficient for simulating the instrument re- 
sponse at a given energy. When simulating the ac- 
tivation of INTEGRAL and SPI by cosmic-ray pro- 
tons, 10 7 protons, requiring a computation time of 
about 3 CPU days (and representing about 21 s of ac- 
tual in-orbit irradiation), provide acceptable statis- 
tics for calculating nuclide production rates, and 
enough statistics in the detector count rate to ob- 



tain a crude spectrum. Detailed studies of prompt 
background line production by cosmic-ray protons 
require multiple simulation runs. 

The MGGPOD web site also provides off-the-shelf 
model spectra for two common radiation fields: 
cosmic-ray pro tons (at solar minumum and maxi- 
mum, based on Moskalcnko et all 120021) . and the dif- 
fus e cosmic X-ray and gamma-ray emission (based 
on lGruber et a l.. 1999). For missions in low-Earth 
orbit albedo radiations (most importantly neutrons 
and gamma rays) and geomagnetically trapped par- 
icles need to be considered in addition. The web site 
provides links to the ESA Space Information System 
(SPENVIS) and the NASA Space Ionizing Radia- 
tion Environments and Shielding Tools (SIREST), 
which can be used to obtain radiation field mod- 
els for missions in a low-Earth environment. Users 
should be aware that the space environment is not 
stable. Detailed simulations require that input par- 
ticle spectra reflect conditions in a given orbit and 
at a given phase of the solar cycle (including the 
polarity of the heliosph ere for cosmic-ray protons, 
iMoskalenko et~afl Effll . Currently, GCALOR does 
not treat hadronic interactions for alpha particles or 
heavier cosmic rays. If diffuse cosmic X-rays are an 
important background component, the uncertainty 
in the simulation result is usually dominated by un- 
certainties in the photoelectric absorption in materi- 
als surrounding the detectors as defined in the mass 
model. 



TGRS Jan. through May 95 Data 
MGGPOD Simulation 




Figure 2. A comparison of the January-May 1995 
TGR S spectrum with a M GGPOD simulation, taken 
from. lU* eidensnointner et all \200A) . The broad fea- 
tures in the data between 210 and 260 keV are elec- 
tronic artefacts 



a strong and time- variable background component. 
Unlike all other components of the radiation envi- 
ronment, the time- variability of the S A A- induced 
backgrounds needs to be accounted for in the simu- 
lat ions. Preliminary mod elling results are presented 
bv lWunderer et al . (2004) at this workshop. 



5. SUMMARY 



4. APPLICATIONS OF MGGPOD 



The MGGPOD Monte Carlo suite has been, and 
is being, applied to model the instrumental back- 
grounds of several gamma-ray missions. Using the 
MGGPOD codes, very good agreement between the 
Monte Carlo results and the actual data of the Tran- 
sient Gamma-Ray Spectr ometer (TGRS) on board 
Wind llOwens et al.lll995ft has been obtained , as dis- 
cussed in Wcidcnspointner ct al. (2003, 2004). As an 
example, in Fig. [5] we show a comparison of MG- 
GPOD simulations with TGRS data. First mod- 
elling results for the SPI Spectro meter on board 
the ESA INTEGRAL observatory IjVedrenne et all 
2003) again yielded good agreement, but also in- 
dicate remaining deficiencies with respect to the 
production and thermalization of secondary neu- 
trons in a massive spacecra ft and/or instrument 
(Wcidcnspointner ct al. , 2003) . Both of these instru- 
ments operate in highly elliptical orbits above the 
Earth's radiation belts. 

Recently, the MGGPOD suite has been applied 
to modelling the instrumental background of the 
Reuven Ramaty High-Energy Solar Spectroscopi c 
Imager (RHESSI, described in ISmith et all 12002^ . 
This is the first time MGGPOD is used for mod- 
elling the instrumental background of an instrument 
in low-Earth orbit. This case is particularly dif- 
ficult because activation during passages through 
the South-Atlantic Anomaly (SAA) gives rise to 



The MGGPOD Monte Carlo suite is an ideal tool for 
supporting the various stages of gamma-ray astron- 
omy missions, ranging from the design, development, 
and performance prediction through calibration and 
response generation to data reduction. The MGG- 
POD software and documentation are publicly avail- 
able for download at CESR. The package has been, 
and is being, successfully applied to several past and 
present gamma-ray missions. 
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